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We show that deep inelastic neutron scattering from hydrogen(or other light nuclei) can be used 
to measure a spectrum of anharmonic contributions to the target atom momentum distribution with 
^\ • high and known accuracy . The method is applied here to determine the momentum distribution of 

On ' the hydrogen in the hydrogen bonded system i('i/C204 (potassium binoxalate), where 13 anharmonic 

On . coefficients are obtained at the 2a to Scr level. The momentum distribution is obtained to an accuracy 

of better than few percent at all significant values of momentum. 
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I. INTRODUCTION 

The measurement of proton momentum distributions by neutron scattering is analogous to the measurement of 



electron momentum distributions by Compton scattering ||^ and measurement of nucleon momenta by Deep Inelastic 
^ ■ Scattering ||^ and is known as Neutron Compton Scattering (NCS) or Deep Inelastic Neutron Scattering (DINS). All 
three techniques rely upon the fact that if the momentum transferred from the incident to target particle is sufficiently 
large, the impulse approximation (lA) can be used to interpret the data. In the lA, momentum and kinetic energy are 
, conserved. From a measurement of the momentum and energy change of the neutron, the momentum of the target 
nucleus before the collision can be determined. 
On DINS measurements have only become practical since the construction of intense accelerator based neutron sources. 

On which have allowed inelastic neutron scattering measurements with energy transfers in the eV region ^J. Energy 
transfers much greater than the maximum vibrational frequency of the target atom are required before the lA can 
be used to reliably determine the momentum distribution . At lower energy transfers the lA is no longer valid and is 
not related in a simple way to the observed scattering intensities. 

There have been a few pioneering studies on anisotropic systems at eV energy transfers but the analysis has 
^ ■ been limited to fitting Gaussians to the observed data, or more generally fitting the data with model containing a 
Q ' few parameters, as was done for measurements on molecular hydrogen We show here that an entire spectrum of 
O anharmonic coefficients can be measured without recourse to any model, in addition to the widths of an anisotropic 
^ ' gaussian, thus describing an arbitrary anisotropic and anharmonic momentum distribution. The possibility of doing 
this for isotropic systems was first suggested by Reiter and Silver, |jl^ That possibility, for more general systems, 
is now a reality. We demonstrate this by measuring the momentum distribution for KHC2O4 where we obtain 14 
anharmonic coefficients whose size varies by nearly two orders of magnitude, with at least 2-3cr confidence levels for 
all but one. The experimental instrument is the EVS spectrometer at ISIS. The work presented here by no means 
represents the limits of resolution of the instrument, but rather the first experiments of this kind. Upgrades are 
planned in the near future that will significantly increase flux and counting efficiency. 



THEORY OF MEASUREMENT 



The theoretical basis of neutron Compton scatte ring is the impulse approximation (lA), which is exact when the 
momentum transfer and energy transfer are infinite |11-13|. The neutron scattering function S{q,Lu), is related to the 
momentum distribution Ti{p) in the impulse approximation limit by the relation 



S{q,uj) = — [ n{p)S{y ~ p.q)dp = —J{q, 



v) 



(1) 



where y=-y (w — ^g), M is the mass of the target particle,q=|(7|, and q—q/q. 

DINS measurements on protons have a particularly simple interpretation, as the interaction of protons with other 
atoms can usually be accurately accounted for |14|-|l6[| in terms of a single particle potential and hence by a proton 
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wave function, [ p^t From elementary quantum mechanics, Ti{p) is related to the Fourier transform of the proton wave 
function via, 

n{p) = \ I / \<i>{r)exp{ip- r)dfp (2) 



(27r)3) 

and a DINS measurement of n(p) can be used to determine the wave function in an analogous way to the determination 
of real space structure from a diffraction pattern. If n(p) is known, and if the proton is in a site with reflection 
symmetry, so that the wave function can be assumed real, then in principle both the proton wave function and the 
exact form of the potential energy well in which the proton sits can be directly reconstructed. |ic| ] With an asymmetric 
site such as potassium binoxalate, the phase information that is lost by taking the absolute value of the momentum 
wavefunction is irrecoverable, and we will not be able to reconstruct the potential directly. The n(p) obtained can, of 
course, be used to check any model potential. 

While the original formulation of the inversion problem [ p^ is complete as it stands, it is useful for the systems we 
will be dealing with to take into account the anisotropy of the system explicitly. The fundamental result that allows 
for a simple inversion of the Radon transform, 3{q, y) to obtain n{p) makes use of a basis of Hermite polynomials and 
spherical harmonics in which the transform is diagonal. That is, a single term in the series for J{q,y) corresponds to 
a single term in the expansion of n(p). 

If we express 3{q,y) in this basis as 

Jiq,y) = — V a,n,l.mH2n+l{y)yira{q) (3) 

7r2 ^—^ 

n.L,m 

then n(p) is given in the related basis of Laguerre polynomials as 

n,t,m 

where p and q are unit vectors. Clearly, since the expansions are complete, a distribution of the form 

2 

— 

(27rtTi)2 

with the ffi significantly different from each other, could be expanded in this form, but even if R(p) were 1, it would 
require a large number of terms in the series. To avoid this, we show that the anisotropy may be taken into account 
by a change of variables, so that the coefficients an^i.m represent genuinely anharmonic contributions. 
Introducing the new variables 

p', = p^/V2a, (6) 
with n{p) defined as in Eq. ^ , defining R' {p')=R{p{p')) and 

_ /2 

n'if) = (7) 
7r2 

we have 

Jiq, y)^ f n\p')5{y - p' .q')dp' (8) 



where q' ^ — qi\f2ai. The right hand side of Eq. (|^) is no longer a Radon transform, since (f is not a unit vector. 
However, defining y' = yl\(i\ we obtain 

^(9, y) = T^ [ n'{p')S{y' - p'.q')dp' -i-J'iq', y') (9) 
W\ J \q'\ 

where J' {q' , y') is the Radon transform of the isotropic(in it's gaussian component) but anharmonic distribution n'{p'). 
If q is specified as a unit vector in the usual spherical coordinates, then 
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\q'\ = V2{{aisin{e)cos{(t)))'^ + {cf2sin{e)sin{(t))f + {as.cos{S)Y)^^ 



(10) 



Our procedure is to expand J'{q',y') in hcrmite polynomials, as in Eq. (|^), and least squares fit the data, ?>{q,uj), 
using Eqs. ([l|,p|,p^ , to obtain the parameters, cri,an,i,m- n'[p') can then be reconstructed as in Eq. (jj), and we thus 
obtain n(p) as in Eq.(|^) with R(p)= R'{p'{p)). That this is a practical procedure will be demonstrated below. 



II. MEASUREMENTS 

The measurements were performed on the electron volt spectrometer eVS [ p^ at the ISIS neutron source. On 
EVS the energy of the scattered neutron is fixed by a resonance filter difference technique [l9| ]. The final neutron 
velocity and energy are related by Ei = where m is the neutron mass. The energy of the incident neutron is 

determined from a measurement of the neutron time of flight via the equation 



t 



(11) 



where t is the measured time of flight , Lq and Li are the lengths of the incident and the scattered flight paths of the 
neutron, and vi are the speeds of the incident and scattered neutrons. Then 



(12) 



and 



(13) 



where 9 is the scattering angle. From these equations lo and q can be determined for a given time of flight t , if the 
instrumental parameters Lq,Li,9 and Ei are known. Hence from the count rate at a given time t, J(q, y) can be 
determined. On eVS the detectors are situated in the horizontal plane and hence q is always horizontal. By orienting 
the sample with a chosen crystal axis vertical, it is possible to measure J(g, y) for g in whichever plane, relative to the 
sample, one chooses. A time of flight scan at a particular angle for a given detector does not correspond, however, to 
a particular direction of q. There is significant curvature of this scan through the proton momentum space since the 
direction of q varies significantly over the data region. Time of flight spectra for eight adjacent detectors at angles 
between 35 and 55 degrees scan through the atomic momentum momentum space of the proton as illustrated in flgure 
1 



SCAN IN Y SPACE, M= 1 .0079E+00 E1= 4,9080E + 03 
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MOMENTUM NORMAL TO BEAM (INV ANG) 
FIG. 1. Scan pattern in momentum space for detectors at a fixed angle as the time of flight is varied 
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A complete scan over the proton momentum space is constructed by combining a number of data sets, taken with 
the sample rotated about the vertical axis by appropriately chosen angles. 

The reported measurements were made using two banks of 8 Li^ doped glass scintillator detectors which were 
symmetrically placed on each side of the incident beam at scattering angles between 35° and 55°. For DINS studies 
of protons it is necessary to site the detectors at forward scattering angles since the hydrogen scattering cross section 
is strongly anisotropic at eV incident energies, with virtually no back scattering. This restriction is a kinematic 
consequence of the closeness of the mass of the neutron and the hydrogen atom and does not apply to heavier atoms. 

The resolution function of the instrument is determined by the uncertainties in the measured values of the time 
of flight t and the distribution of L(),Li,6 and Ei values allowed by the instrumental geometry and analyser foil 
resolution. Uncertainties in Lq arise primarily from the finite depth of the neutron moderator, those in Li and 
6 from the finite sample and detector sizes and those in t from jitter in the detector electronics. All resolution 
components can be determined by calibration measurements and all except the energy component can be approximated 
by Gaussians, without significant error. A 0.015 mm thick gold foil provided a Lorentzian energy resolution function at 
El = 4908rney , with a peak HWHM of 136 meV. The Gaussian and Lorentzian resolution components in momentum 
space y ,are listed in table 1 for two angles representative of the range of angles employed. The resolution is dominated 
by the energy component which varies strongly with scattering angle. The second most important contribution comes 
from the angular resolution of the spectrometer and is independent of angle. The momentum and energy transfers at 
the centre of the hydrogen response peak are also listed for the different angles. 

Table 1. The resolution widths are the Lorentzian HWHM for (RL) and the Gaussian standard deviation for other 
parameters (RG) . The momentum q and energy transfer ( at the scattering angles 35° and 55° are also given.) 



Angk" 






q ^-') 


^' (eV) 


35*^ 


0.61 


1.08 


34.1 


2.41 


55^ 


0.55 


0.55 


48.8 


4.92 



The raw data contains signals from all the atoms in the scattering sample and from the cryostat background. Fortu- 
nately the energy transfer to hydrogen is much greater than that to other atomic masses and the proton signal is well 
separated from that due to other masses. The contribution from all components other than hydrogen is subtracted 
by fitting a sum of Gaussians convoluted with the instrument resolution function to the data and subtracting off the 
fitted contribution to other peaks. There is also a small contribution to the data from a second gold resonance at 60 
eV which can be seen at 100 /xsec and this is also fitted and subtracted from the data. 

Time of Flight Data 

3x10"^ 

2x10"^ 
c 







200 400 
Time of Flight (fi sec) 

FIG. 2. The sum of data from 8 detectors at scattering angles between 38° and 55° is shown as the dotted line. The data 
after subtraction of the contribution from atoms with higher masses and the 60 eV resonance data is shown as the full line. The 
total data set for a single plane consisted of 36 such spectra. The FWHM of the Lorentzian resolution function is 6 microsec. 

The data for each scan was converted into a distribution in the momentum space of the crystal as described above. 
Complete coverage of the plane was achieved by combining six runs which were taken at steps of 23° sample rotation. 
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A contour plot of the data derived from the six runs is shown in figure 3. This was produced by binning the counts 
from different detectors and sample orientations in the appropriate pixel of the momentum space 
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FIG. 3. Contour plot of data in one plane of potassium binoxalate. The apparent double peak near the origin is an artifact 
of the method of plotting the data 

The data has been corrected for sample attenuation, but still contains errors due to small deviations from the 
Impulse Approximation which are present at the finite momentum transfers of the measurement. These tend to 
introduce small asymmetries into the data set at the 1-5% level, thereby removing the exact inversion symmetry of 
the Compton profile. It has been shown by Sears that most of these effects are removed by symmetrisation of the 
data about the origin. This procedure cannot be followed for our data sets as the scan in a single detector is curved 
and different points in the crystal plane which are related by inversion symmetry may have been measured in different 
detectors, with different experimental resolution. However by fitting to a, 3{q,y with inversion symmetry, as discussed 
below, we automatically include a correction for deviations from the impulse approximation. Any asymmetries are 
ignored by the fit and should not affect the values of the fitted parameters. 

III. FITTING PROCEDURES 

Eqs.(^j^) hold quite generally for Radon transform pairs, but physical requirements in the present context restrict 
the allowed coefficients. Since J((f, y) is an even function of y, 1 is restricted to even values, and since J{q, y)is real, 
the an^i^m for ±m must be equal. In fact, if there are residual final state effects in the data, the data will not be 
symmetric. Restricting the coefficients in this way will therefore eliminate these residual effects. Lifting the restrictions 
and fitting the data allows one to measure the extent of these effects. For potassium binoxalate, data was taken for 
three perpendicular planes oriented parallel to the crystal axes The procedure followed was to perform a simultaneous 
fit to the 6x16x3=288 separate time of flight spectra to an expansion of the form given in Eq.(^, convoluted with the 
instrument resolution function 

The actual fitting procedure requires that three of the parameters Gi, an,i,m be fixed, as they are not all independent. 
That is, ff 0116 Vcirics tliG (Ji cibitrarily, there will alwRys be cl set of (in,i,m 

that will fit the data, the particular values 
that fit depending on the choice of the Ui. One could obtain fixed values for the ai by first fitting the data to a gaussian, 
i.e. R'^j?) = 1, and then varying only the an,i,m, in which case all such terms would in principle be needed. We find 
in practice, that better results are obtained, i.e. fewer anharmonic coefficients needed and more rapid convergence of 
the series for R'ijf), when the ai are allowed to vary but the first three anharmonic coefficients, ai_o,0) «o,2,Oj ^0,2,2 are 
set to zero. In fact, these coefficients are not really anharmonic coefficients at all. They could always be eliminated 
by a shift of the ai and an adjustment of higher order coefficients. This procedure thus has the virtue as well of 
producing only genuinely anharmonic corrections to a gaussian fit. 

The datasets as they are presently obtained in the EVS spectrometer at ISIS, are obtained one plane at a time(see 
discussion of experimental apparatus). That is q varies within a plane, and a range of y values is taken such that 
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J((f, y)is negligible outside this range. There is a very high density of points, which for the purposes of the present 
discussion we can take to be continuous. The question then arises, 'How many planes of data are needed to determine 
a specified number of coefficients, and at what angles to each other should they be?' We can see the point of the 
question by considering first a single plane and looking at the fit to the leading anharmonic coefficient. From Eq.(|^) 
we see that there are six independent coefficients multiplying H^, i.e. the coefficients of (1001^^207^22,^401^427^14). 
Let us say that our coordinate system is chosen so that the plane measured is treated as an xz plane. Then in terms of 
the variable 9 there are only three independent fourier coefficients that can be present in the data for the coefficients 
of H4, i.e. the coefficients of (f , cos{29), cos{A9) for instance. Therefore three of the coefficients are not independent. 
The complete set of coefficients cannot be determined by the data. This is of course due to the fact that there is 
no information as to the behaviour of J{q, y) for nonzero azimuthal angles (p in the data, so we shouldn't expect the 
fitting procedure to provide it. The data provides a complete description of J(g, y) only if this function is rotationally 
invariant about the z axis. Actually, what is required is only that J'{q',y') be rotationally invariant, since the fit is 
done in the primed coordinate system. If this is the case, then all coefficients of Yim with non-zero values of m must 
be zero. We see that there are only three remaining possibly non zero coefficents, which can all be determined. For 
higher order terms as well, keeping only the coefficients with m=0 provides all the independent terms needed to fit 
the data, and the resulting fit, of course, is rotationally symmetric about the z axis. 

If the data is not known to be rotationally symmetric, additional planes of data must be taken to determine even 
these lowest order coefficients. In general, whenever we take another plane of data, we might expect to obtain 3 more 
independent measurements of the coefficients of i?4, 4 independent measurements of the coefficients of H^,^ and in 
general, k+l measurements of the coefficients of i?2fe- (k+1 being the number of independent fourier components in 
the data for that value of k). Since the number of an.i,m that are to be determined for 2n-|-l=2k is (k+2)(k-|-l)/2, it 
appears that (k+2)/2 planes are needed to measure all coefficients up to i?2fc- The angles between the planes must 
be chosen, however, so that the measurements are really independent. For instance, if k=2, it would appear that two 
planes would suffice, but if they are chosen as the xz and yz planes, they do not provide independent measurements of 
the coefficients. This may be seen by observing that the sum of the data from the two planes gives three independent 
fourier coefficients to determine four independent a„^;^m, the coefficients of (5^007^207^07^4)7 since the coefficients 
of Y22 and I42 cannot affect this sum. The difference of the data on the two planes gives three equations for the 
two coefficients of I22 and 1^42. A better choice for the planes would be (/) = and — 7r/4, which would allow the 
determination of all the coefficients. If there is some symmetry in the problem, one may be able to use perpendicular 
planes if the symmetry axis is chosen appropriately with respect to the common axis of the two scattering planes. For 
instance, if there is tetragonal symmetry present, and the symmetry axis is chosen perpendicular to the common axis, 
one can obtain all the allowed coefficients up to k=4. One can show that three perpendicular planes do in fact suffice 
to determine the coefficients up to k=4 in the general case, without any symmetry to reduce the number of allowed 
coefficients. The question of whether this is an optimum configuration of planes or not, we will leave to another time. 
Including the three cr^, a three plane measurement allows 34 coefficients to be measured. 



IV. MEASUREMENT ERRORS 



The uncertainty in the measurement of n(p) at some point p is due to the uncertainty in the measured coefficients. 
Denoting an arbitrary coefficient by p^, we have 

5n{p) = Y: ^-5,, (14) 

The fitting program, after a minimum is obtained with some set of coefficients, calculates the correlation matrix 
< SpiSpj > by varying the coefficients slightly and calculating the curvature of the chi-square of the fit. |^l|. Hence, 
the variance in the momentum distribution is 

There are of course, potential systematic errors that could enter the measurement, such as multiple scattering effects, 
or an error in determining the resolution function. The former must be handled with good experimental design, and 
if small, can be corrected for. We have done measurements on samples whose thickness differed by a factor of two, 
with no significant differences in the observed scattering. Multiple scattering effects would also lead to asymmetries 
in 3{q,y) that are not observed. The resolution function has been studied extensively, and is believed to be known 
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accurately.In the present measurement, the resolution width is between 15 and 25% of the width of the distribution 
measured. 

A further source of error not contained in the estimate in Eq. ( [Tsl ) is the truncation of the series used to fit the data 
to include only terms up to 2n+l=8. To get an idea of the seriousness of this, and to test the fitting procedures and 
software, we have generated synthetic data from a known momentum distribution that corresponds to an asymmetric 
double well, and convolved it with the instrumental resolution function. The data was then analyzed by the means 
described above, and the extracted n(p) compared with the input. The input n(p) corresponded to a spatial wave 
function consisting of two displaced gaussians with the same variance and a relative weight r=.5. The explicit form is 



(1 + H + 2rcos(2p.a)) -PT e 
n(Px,Pv,Pz) = ^ I I -T (16) 

(l + r2 + 2re"^) t 

where Ux = 4.6, CTj, — 4.0, ct^ = 6.0 and a=.15 in units of inverse angstroms and angstroms, respectively. This 
form is rather similar to the actual form of the data we will analyze. The coordinate system z axis was chosen to 
be identical with the crystal z axis. The comparison is shown in fig 4. The extracted n(p) is plotted with a dotted 
line, the input n(p) above with a dashed line. As may be seen from the figure, there is essentially no error due to the 
truncation of the the expansion. Of course, an input with more variation might require including higher terms in the 
expansion, and hence taking more planes of data. The actual data we have obtained for potassium binoxalate has less 
variation than the simulation and the truncation error will be negligible. In the other directions, the fit is rigorously 
gaussian, since all coefficients with m 7^ were zero. We note that the measured parameter was 4.47, not 6.0. It is 
simply a fitting parameter. The overall momentum distribution has physical significance, the individual parameters 
may not. In fact, we have gotten the same degree of fit with the crystal z axis aligned along the coordinate y axis. In 
this case, all the coefficients are non-zero, but the resultant n(p) is identical to the one displayed above. 



asymmetric double we 



20 



25 



FIG. 4. Comparison of input n(p), given by Eq.(16)with reconstruction of n(p ) using fitting procedure. The z coordinate 
axis is chosen to be the double well axis. 



A final source of systematic error involves the possibility of finding a false minimum with the fitting procedure. 
With such a large number of parameters, there is the possibility that the program will home in on a local minimum and 
miss the true minimum. We are using straightforward gradient methods for the search, which would have difficulty 
with a very rugged chi-squared landscape. We do not appear to have such a landscape for these problems, and 
certainly not for the large parameters, such as the at, or the largest anharmonic coefficients, whose values appear to 
be quite robust with regards to different paths to the minimum. One can check for this problem by orienting the 
fitting coordinate system differently with respect to the crystal axes, as was done for the test data above. In this 
case, all the coefficients will be different, but the final n(p) must be the same. There is also the fact that n(p) must 
be positive, and a spurious fit that leads to significant negative values can be rejected. We can also eliminate some 
parameters which are consistently much smaller than their variance, and which can be set to zero without affecting 
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significantly the chi-square, thus reducing the dimensionahty of the space. We beheve there are no problems of this 
sort with the fits we will present. 

V. RESULTS FOR POTASSIUM BINOXALATE 

Potassium Binoxalate, KHC2O4, is a hydrogen bonded system in which the hydrogen sits in an asymmetric position 
between two oxygen atoms. |2^.The crystal is monoclinic. We will choose a primitive cell for which the bond axis, 
that is the line joining the two oxygen atoms, is essentially aligned with the c axis of the crystal. We choose this 
axis for the z axis of our coordinate system. Three planes of data were taken at right angles to each other, with 69,632 
data points in all. One of these is the be plane, where the b axis is the unique axis, and the other two are the a*b and 
a*c plane, where the a* axis is perpendicular to the be plane. These were fit with the methods described above to 
give the momentum distribution. The measurements were done at 10 deg K and hence there are no significant finite 
temperature corrections to the ground state momentum distribution due to excited states. 

Note that the coordinate system used to describe the momentum distribution need not have the symmetry of the 
crystal. The local symmetry of the site is only a two-fold rotation. We show the values of the fitted coefficients in 
tables 2 and 3, along with the rms uncertainty in their values. This is included only to give some sense of the individual 
parameter. The error in n(p)is given by Eq.(Q) and includes the effect of the correlations between coefficients, whereas 
the figures cited in the table are only the diagonal correlation coefficients. It can be seen, that many of the measured 
coefficients have been determined at the 2-3 a level of confidence, some at higher levels, and only one at a la level. 
Coefficients that are set to zero in the fitting procedure were found to have values smaller than their variance by 
at least a factor of 2, and setting them to zero did not significantly change the minimum value of chi-square. The 
goodness of fit to a sample of the data is shown in Fig. 5. 

Constant Angle Plot 




3^ | )O0C00C00C i 



-20 20 

R 

FIG. 5. Data fitted by method described above. The data is a composite of data points in a lOdeg wedge about the z axis. 

We have compared the fitted prediction, convolved with the instrumental resolution, and the data, for the cumulative 
sum of data in a 10 degree wedge along the hydrogen bond direction. The measured momentum distribution along 
the axes of the measurement planes, are shown in Figs. 6-9, and contour plots along coordinate planes in Figs. 10-13. 



8 



potassium binoxalate 




5 10 15 20 



z axis 

FIG. 6. Momentum distribution for potassium binoxalate along the axis shown. The momentum is in units of A~l, n(p is 
in arbitrary units. The errors are calculated as in Eq.(15) for the parameters that are significant in Tables 2 and 3. The lower 
curve is the anharmonic component. 



potassium binoxalate 




FIG. 7. Momentum distribution for potassium binoxalate along the axis shown. The momentum is in units of A~l, n{p) is 
in arbitrary units. The errors are calculated as in Eq.(15) for the parameters that are significant in Tables 2 and 3. The lower 
curve is the anharmonic component. 
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potassium binoxalate 



z 




y axis 

FIG. 8. Momentum distribution for potassium binoxalate along the axis shown. The momentum is in units of A~^, n(p) is 
in arbitrary units. The errors are calculated as in Eq.(15) for the parameters that are significant in Tables 2 and 3 
. The lower curve is the anharmonic contribution. 
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FIG. 9. Momentum distribution for potassium binoxalate along the axes shown. The momentum is in units of n(p) is 
in arbitrary units. 
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oxalic acid 




z axis 

FIG. 10. Momentum distribution for potassium binoxalate along the axes shown. The momentum is in units of 



oxalic acid 




X axis 

FIG. 11. Momentum distribution for potassium binoxalate along the axes shown. The momentum is in units of . 
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potassium binoxalate anharmonicity 




FIG. 12. Anharmonic contribution to the momentum distribution for potassium binoxalate along the axes shown. The 
momentum is in units of A'^. 



potassium binoxalate anharmonicity 




FIG. 13. Anharmonic contribution to the momentum distribution for potassium binoxalate along the axes shown. The 
momentum is in units of A~^. 
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The coefficients that were measured are given in the tables below. 

Table 2. The harmonic fitting coefficients and the variances in their values as measured 



Harmonic Coefficients 


i 




5(Ji 


X 


4.509 


.0235 


y 


4.869 


.0254 


z 


5.351 


.0449 



Table 3. The anharmonic fitting coefficients and their variances as measured. 



Anharmonic Coefficients 


n 


1 


m 


/ 




2 








0.0000 


.0000 


1 


2 





-0.1443 


0.0377 


1 


2 


2 


0.0000 


0.0000 





4 





0.0510 


0.0086 





4 


2 


-0.0564 


0.0095 





4 


4 


0.0000 


0.0000 


3 








-0.1356 


0.0098 


2 


2 





-0.0296 


0.0105 


2 


2 


2 


0.0000 


0.0000 


1 


4 





-0.000 


0.0000 


1 

j_ 


A 


9 


n nnnn 

u.uuuu 


n nnnn 

u.uuuu 


1 


4 


4 


0.0000 


0.0000 





6 





-0.0184 


0.0065 





6 


2 


-0.0127 


0.0028 





6 


4 


-0.0388 


0.0129 





6 


6 


0.0000 


0.0000 


4 








-0.0290 


0.0029 


3 


2 





0.0000 


0.0000 


3 


2 


2 


0.0000 


0.0000 


2 


4 





-0.0016 


0.0013 


2 


4 


2 


0.0069 


0.0013 


2 


4 


4 


0.0000 


0.0000 


1 


6 





-0.0049 


0.0021 


1 


6 


2 


0.0000 


0.0000 


1 


6 


4 


-0.0098 


0.0036 


1 


6 


6 


0.0000 


0.0000 





8 





0.0000 


0.0000 





8 


2 


0.0000 


0.0000 





8 


4 


0.0000 


0.0000 





8 


6 


0.0000 


0.0000 





8 


8 


-0.0075 


0.0028 



Note that there are 14 anharmonic coefficients that are measureable, 13 of which are at the 2u level at least, and 
that the harmonic coefficients are measured to better than 1%. 



VI. CONCLUSION 

We have shown how the method of analysis of DINS data suggested in can be extended to anisotropic momentum 
distributions, and have applied this method to an analysis of the hydrogen bond in potassium binoxalate. The results 
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demonstrate that DINS, as it is implemented now at ISIS, is capable of detailed, model independent, measurement of 
the momentum distribution for hydrogen, and by inference, other light atoms. These measurements required about 
four days of beam time, and are the first such measurements to be analyzed in this way. The count rates, resolution 
and detetector efficiency can be, and are scheduled to be, significantly improved. The data can be analyzed in less than 
a day. The DINS technique thus provides a practical means of accessing precise information about the anharmonicity 
of local potentials and can provide a check of any theoretical calculation of these potentials at a level of accuracy and 
detail that has not been possible previously. 
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APPENDIX A: SEPERABLE DISTRIBUTIONS 



The distribution used for the simulations described here is a special case of a general separable distribution 

nip}^Y[n,{p,) (Al) 



If we represent 



make use of the fact that 



6 J) 

= 77 TT y] ai.nHn(— ^) (A2) 



e*«^-^'iJ„(a;)da; = 7r3e^(ig)", (A3) 



and represent the delta function in the definition of the radon transform by its fourier transform, it is straightforward 
to show that the Radon transform of n(p) in the isotropic coordinate system defined in the text is given by 

J'{q\y') = ^E(n««^".(4)"0^.H+«.+n3(y') (A4) 

For the case we have used as a simulation, in which the distribution is anharmonic only in the z direction, the result 
simplifies to 

J'{q',y') = ^5](a3,„(co5(0')")^«(y') (A5) 

In this case, the coefficients in the expansion of J'{q', y') in the form given by Eq.(^re not all independent. There is 
only one independent coefficient in the set of a„_;_o for each value of 2n+l. 
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